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Abstract. 

C\| | We confront Cardassian models with recent observational data. These models can 

be viewed either as purely phenomenological modifications of the Friedmann equation, 
or as arising from cosmic fluids with non-standard properties. In the first case, we find 
OO ' that the models are consistent with the data for a wide range of parameters but no 

CD ■ significant preference over the ACDM model is found. In the latter case we find that 

the fact that the sound speed is non-zero in these models makes them inconsistent with 
the galaxy power spectrum from the Sloan Digital Sky Survey. 
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1. Introduction 

o 

The Universe is expanding at an accelerating rate. We have direct evidence for 
this from the observed magnitude-redshift relationship of supernovae type la (SNIa) 
PU |2] • Cosmic acceleration requires a contribution to the energy density with negative 
^ ■ pressure, the simplest possibility being a cosmo logical constant. Independent evidence 
for a non-standard contribution to the energy budget of the universe comes from e.g. 
the combination of the power spectrum of the cosmic microwave background (CMB) 
temperature anisotropies and large-scale structure (LSS): on the one hand, the position 
of the first peak in the CMB is consistent with the universe having zero spatial curvature, 
which means that the energy density is equal to the critical density. On the other 
hand, observations of large-scale structure, e.g. the clustering of galaxies, show that the 
contribution of standard sources of energy density, whether luminous or dark, is only a 
fraction of the critical density. Thus, an extra, unknown component is needed to explain 
the observations Ej • 

The model resulting from reintroducing Einstein's cosmological constant, A, in a 
universe with baryons, radiation, and cold dark matter (CDM) is consistent with all 
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large-scale cosmological observations like the anisotropies in the CMB radiation and 
the power spectrum of galaxies [5J. However, the modern interpretation of A is that it 
represents the energy of the vacuum, and from the viewpoint of quantum field theory, 
the tiny value inferred for it from observations is highly unnatural [Hj. Faced with 
this problem, the popular choice is to set A to zero and invoke a new component to 
explain the acceleration. This does not solve the problem of the smallness of A; if it 
is indeed equal to zero, one still needs to understand the physical mechanism behind 
this. However, one may hope that A = may be easier to explain than a small, but 
nonzero A. The question is then what the unknown component driving the accelerated 
phase of expansion is. One needs to introduce a component with negative pressure, 
and this can be done e.g. by invoking a slowly evolving scalar field [3 El Ej , or a 
negative-pressure fluid, e.g. a Chaplygin gas [HUE] A negative-pressure fluid, however, 
can be problematic due to its fluctuations. Fluctuations of the unknown fluid can grow 
very rapidly and hence LSS surveys can place strict constraints on such models (see e.g. 
[l2l[T3]). In order to circumvent this, it is sometimes assumed that the new fluid does 
not fluctuate on the scales of interest and its existence is visible only through modified 
background evolution. 

As an alternative to fluid models of dark energy, modifications of the Friedmann 
equation have been considered by several authors [HlHIIllinilUllHllIHllini- m an earlier 
paper [21], we investigated a model E] where the left-hand side of the Friedmann 
equation is modified. In the present paper, we will look at the so-called Cardassian 
models ^11201; where the right-hand side of the Friedmann equation is modified. 

The structure of this paper is as follows. In section 2 we summarize the basic 
properties of the Cardassian models. In section 3 we consider fits to SNIa and WMAP 
data. In section 4 we extend the analysis to include probes of the large-scale distribution 
of matter, and in section 5 we summarize our constraints. Finally, section 6 contains 
our conclusions. 

2. Cardassian models 

Cardassian models were introduced in (TH] as a possible alternative to explain the 
acceleration of the universe by a model that has no energy components in addition 
to ordinary matter. Necessarily, the Friedmann equation must be modified: 



where p m is the energy density of matter (baryonic and dark), p r is the energy density of 
radiation radiation and n, A are parameters of the model. The requirement that universe 
undergoes acceleration at late times sets n < 2/3. Phenomenologically there is no lower 
limit but conventionally only models with n > are considered. The appearance of the 
extra term on the right hand side was originally motivated by brane cosmology [22J. 
This approach was, however criticized in (22] where it was shown that in the scenario of 
[22]; the weak energy condition was violated on the brane. We will ignore the radiation 
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term in the following discussion for simplicity but keep it in the CMB calculations as 
there it will have an effect. For SNIa constraints, the radiation contribution is negligible. 

More recently, the Cardassian scenario has been generalized (see e.g. [20]) and its 
connection to fundamental physics has become less apparent. Instead, one can view the 
Cardassian models as a class of phenomenological models that attempt to explain the 
accelerated expansion of the universe without explicitly resorting to adding new cosmic 
fluids, such as a cosmological constant, to the Friedmann equation. 

A most general formulation of the Cardassian framework is to consider Friedmann 
equations that depend only on the matter density in the universe, thus 

H 2 = yffW, (2) 

where p m is the energy density of matter (CDM+baryons). Nucleosynthesis dictates that 
at early times, the evolution of the universe must be very close to the standard H 2 ~ p, 
hence g{p m ) ~ Pm for large p m . In principle, one can consider arbitrary functions g or 
expansions such as 

oo 

g(f>m) = Pm + J^CiPm > ( 3 ) 

i=0 

where q are parameters. The original Cardassian scenario then corresponds to 
considering just the leading term in the expansion. 

A commonly considered model within the Cardassian framework is the Modified 
Polytropic Cardassian (MPC) model [20], which extends the original Cardassian 
proposal by introducing an additional parameter. In the MPC model, 

SCO = 9 MPC ( P J = Pm (l + , (4) 

v Pc 7 

where q is the new parameter. The original scenario is reached in the special case q = 1. 
Specializing to a flat universe, and defining Q m = 87rGp m /3ifo (where p m is the present 
day value), the MPC model can be written as 

' ' 2 = n m a- 3 (l + Oca 3 *-"') (5) 



>H 

where VLq = (fi m 9 — 1)- Similarly, the generalized model, Eq. |2l can be written as 



/(«), (6) 



with f(a) in the MPC case given by Eq. (jSJ). 

The main observational constraints to any non-standard scenario are SNIa, CMB 
and LSS. Which probe is most effective depends on the details of the considered scenario. 
If the model in question contains new, fluctuating fluids, the gravitational clustering can 
be radically modified [1] and can lead to very tight constraints. If, on the other hand, 
only the background evolution is modified via the Friedmann equation, the growth of 
large scale structure, while modified, may not differ too radically |2"H 125]. 

The Cardassian scenario has two possible physical interpretations: the first and 
maybe most common one is to consider the modified Friedmann equation arising from 
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modified gravity/extra dimensional constructions, leading to effective terms on the RHS 
of the equation. Another alternative [201; the so-called fluid interpretation, is to consider 
the modified Friedmann equation as a result of having an additional exotic fluid in the 
universe. 

Here, we consider both possibilities using data from SNIa, CMB and LSS 
observations to see how well the Cardassian scenarios are constrained by the combination 
of current data. First we study the effective model with no fluctuations in the Cardassian 
fluid and later consider the effects of the fluctuations. 



3. Fits to CMB and SNIa data 

3.1. Modified background evolution 

When fitting the CMB data with the Cardassian models, one should in principle start 
from the full underlying theory. This can be problematic. If e.g. the Cardassian models 
arise from a scenario with extra dimensions, studying the evolution of perturbations is 
difficult because of the bulk-brane interactions and hence we will in the following use 
a simplified approach where we solve the standard 4-d perturbed Boltzmann and fluid 
equations, but with the background evolution given by Eq. (JTJ) (also see discussion in 

For convenience, we will parameterize the effect of the extra term in the Friedmann 
equation by a dark energy component with an effective equation of state w(a). As long 
as the extra component does not fluctuate, i.e. it only has an effect on the background 
evolution, such a parameterization is equivalent to modifying the Friedmann equation. 

For a general Friedmann equation, Eq. (|2j), we can write the expansion of the 
universe in terms of an effective dark energy fluid. A dark energy component with 
equation of state w(a) has a density which varies with the scale factor a according to 

f f a da' } 

Pd(o) = PoDexp <! —3 y —[l + w(a'))j, (7) 

so for a flat universe, — 1 — the standard Friedmann equation is 

(i) 2 ^ + ( i - n »)^{- 3 r^ ii+t " (a,)| }- <8) 

Using this along with Eq. (|B|). it straightforward to see that the effective equation 
of state is (for f\a) ^ Q m a~ 3 ) 

a f'(a) + 3fi m a" 4 

w (a) = -!--— — — — . (9) 

S j (a) — \l m a 6 

For the special case f(a) = Q m a~ 3 + (1 — Q m ) corresponding to the ACDM model, 
w(a) = —1 as expected. In the case of the MPC model, the equation of state of the 
effective dark energy fluid is 

a 3(l-«k fic + fica 3 9 (l-nA 1/9 1 

w(a) = (l-n) ^ ' . (10) 

i + n c a^- n )A 
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Figure 1. The equation of state w(a) as a function of log(a), for f2 m = 0.3. 
Left panel: n — 0.2, q = 1 (top), 3, 5, 7, 9 (bottom), right panel: q = 3, n = 
0.0 (bottom), 0.1, 0.6 (top). 



At early times, a — > 0, the equation of state tends to a constant: w — > q(n — 1). Also 
note that for the special case q = 1, the equation of state is constant w = n — 1. For 
q ^ 1, w evolves in such a way that w is more negative in the past and increases at late 
times. The equation of state of the effective fluid in the MPC model is shown in Fig. ^ 
for different values of the parameters. 

3.2. SNIa and CMB data 

For the supernova fit we use the sample of 194 SNIa presented in [21]. We fit three 
parameters to to the data: f2 m , n and q and values considered are < Q m < 1, 
< n < 2/3 and q > 1. The Hubble parameter h is also involved in the fits, but is of 
little interest here, and we choose to marginalize over it (by maximizing the likelihood). 

For fitting the CMB TT power spectrum j2H] we use the likelihood code provided by 
the WMAP team [29j||. The CMB power spectra is computed by using CMBFAST code 
[SHj, version 4.5. l^f. For each point in the parameter space (f2 m , n, q, h), we calculate the 
CMB TT power spectrum keeping the amplitude of the fluctuations a free parameter 
by finding the best fit amplitude for each set of parameters. In other words, we only 
fit the shape and not the amplitude of the power spectrum. In calculating the CMB 
power spectrum we use fib = 0.044 (so that we vary the density of cold dark matter, Q c ) 
and ignore reionization effects. We have checked that setting the reionization optical 
depth to the WMAP value makes very little difference in the resulting contours. Fixing 
fib/fim instead of does change the contours somewhat. For example in the origincal 
Cardassian model the CMB contours expand to include higher values of n but the Q m 
range within the contours stays the same. 

Parameters fi m 6 [0, 1] and n G [0, 2/3] have uniform priors and are chosen to cover 
the most interesting range of parameters. The typical grid size used in the parameter 

|| http:// lambda.gsfc.nasa.gov/ 
http://www.cmbfast.org/ 
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scan was typically (h,Q m ,n,q) ~ 5 x 30 3 (larger grids were also used to scan larger 
parameter spaces than what is shown here). For the Hubble parameter h we use a 
Gaussian prior based on the HST Key Project value h = 0.72 ±0.08 [3T] and marginalize 
over h in making all the plots unless otherwise stated (we find that the choice of a flat 
or a Gaussian prior makes little difference). 

3.3. The original Cardassian model 

The original Cardassian model ^H] corresponds to the choice q = 1, leaving only three 
free parameters (fl m , h, n). This model has been extensively studied by using SNIa data 
[H2J EHJ EU HHH HUH HH] ■ The overall picture is that supernovae alone cannot constrain 
the parameter space particularly well and leave quite a large degeneracy in the (fl m ,n) 
plane. CMB data has been utilized in [SHEHj, but only using the power spectrum peak 
positions. Finally, in [39J, the original Cardassian scenario has been constrained by 
Sunyaev-Zeldovich/X-ray data. 

Here, we use the full shape of the WMAP power spectrum along with SNIa 
observations as well as consider the more general MPC model. In Fig. 121 we show 
results for the original Cardassian model from fitting to the SNIa and WMAP data. We 
see that compared to the SNIa data, the WMAP data constrains the model much more 
effectively. The combined contours demonstrate how the original Cardassian model is 
quite severely constrained by the current data. The original Cardassian model seemingly 
does not offer any advantage over the concordance ACDM model. 

3.4. The MPC model 

In the case for the more general MPC model, the situation is however somewhat different. 
The combined fit to CMB+SNIa in the model is shown in Fig. El and the different 
sub figures show the confidence contours after marginalizing over one of the remaining 
parameters f2 m , n, q. From the figures we see that the data prefers fl m ~ 0.3, which is 
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Figure 3. Joined 68, 95 and 99% confidence contours for Q m and n (left), f2 m and 
q (middle), and q and n (right), from fitting to SNIa and CMB data. In each figure, 
we have marginalized the Hubble parameter with a HST prior and the remaining 
parameters with a uniform prior. 



expected in light of the tight constraints coming from the WMAP data. Marginalizing 
over Q m , we see that there is quite a large degeneracy in the (q, n) plane. Comparing to 
previous works [20, 40J however, we see that utilizing the full CMB spectrum along with 
SNIa data tightens the constraints significantly. Even though in the full MPC model the 
data prefers a non-standard model, there is no significant preference over the ACDM 
model (q = 1, n = 0). 

In a recent work |H] , the CMB TT power spectrum in the Cardassian scenario was 
also considered and it was shown that the large scale part of the spectrum is modified, 
on par with what we find here and what we expect based from analytical approximations 

4. The Fluid Interpretation 

The Cardassian models can have different interpretations depending on the properties 
of the fluid. One can either consider that the modified Friedmann equation arises from 
modified gravity, e.g., in extra dimensional theories, or that gravity is standard but 
the fluid has unusual properties. The fluid interpretation was introduced in [20] as a 
possible scenario responsible for the modified Friedmann equation. 

The basic idea of the fluid interpretation of the Cardassian models is that gravity 
is four dimensional but the matter fluid has some pressure which leads to a Friedmann 
equation that is effectively modified. The starting point is the set of ordinary Einstein's 
equations, 

R^-^un = 8nGT^, (11) 

where the energy momentum tensor has the ordinary perfect fluid form T^ u = 
diag(p, — p, —p, —p). The Friedmann equation is of the normal form H 2 = (8nG)p/3. 
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The crucial point is that now p ^ p m but p = p(p m ) instead. Energy density must be 
conserved, i.e., Tj£ = 0, which leads to the usual continuity equation for the fluid p: 

p + 3H(p + p) = 0. (12) 

In order to have the modified Friedmann equation, Eq. (J2J), we must make the 
identification p = g(p m ). Using this in the continuity equation along with the fact 
that the matter fluid has no pressure, gives an expression for the effective pressure of 
the fluid 

P = Pm9'(pm) -g(Pm), (13) 

where ' = d/dp m . It is then straightforward to calculate the sound speed (in the 
adiabatic approximation) and the effective equation of state of the fluid: 

Cl = Pmy (14) 
U = Pm— - 1- (15) 

9 

Physically these properties require that the matter fluid must have some interactions 
that lead to the pressure term of the form Eq. (jT3Jl (see (201 for a discussion). 

Here we only consider the effects on the growth of large scale structure as it is a 
good discriminator for models with a new fluctuating component. The CMB spectrum 
was studied in [41], where they find that the large scale part of the spectrum can be 
modified radically in the MPC scenario with fluctuations. 

4-1. Growth of perturbations: general case 

Since the effective fluid has non-zero pressure, perturbation growth can be dramatically 
different compared to CDM. 

For a single fluid, the fractional perturbation of the effective fluid 5 = (Sp)/p, in 
momentum space obeys [IB] (again in the adiabatic approximation, for modes for which 
k^H) 

5 k + H(2 - 3(2uj - c*))4 - *-H 2 (l - Qc 2 s + 8uj - 3u 2 - = 0. (16) 

The physically interesting quantity is the perturbation of the matter fluid, 8 m , which 
is related to 5 by 

5 = ^ = p m -5 m = (l + ^)5 m , (17) 
P 9 

from which it is clear that the two perturbations agree only when the fluid has no 
pressure (w = or g(p m ) = p m ). 

From Eq. one finds that the equation describing the evolution of the matter 
perturbation is 

5 m ,k + H(2 - 3c 2 J m , k + { - hi 2 (l - 2c 2 s - 3c 4 s + 9(c 2 s - u) 2 ) 

+^ + rz-K* = ' ( 18 ) 

a 2 1 + uj ) 
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where we have made use of the identity 

u = -3H(l + u)(<%-u). (19) 
With c 2 = 0, lj = (or g(p m ) = p m ) this reduces to the standard CDM expression, 

5 m , k + 2H5 m , k - ^H 2 5 m , k = 0. (20) 
4-2. Growth of perturbations: MPC model 

Specializing to the MPC model where g is given by Eq. (J3J, we find the following 
expressions: 

u = (n - 1) 1 „ (21) 

2 . Y-, , ( ra - !)? 



c f = ^ 1 + 1 ' 'I » , (22) 

where have used the previously defined quantity Qc = (^m 9 — !)■ From the expressions 
we see that in the fluid interpretation, the ACDM model is reached in the limit 
n — 0, g=l. 

We have solved the equation describing the growth of the matter perturbation, 
Eq. (|18jl . numerically in the MPC model. Following |13j . we use CMBFast to calculate 
the matter transfer function from an initial scale free spectrum in the ACDM model 
at z = 200 and then evolve the perturbation according to Eq. (fTSJ). The resulting 
matter power spectrum is shown in Fig. E]for different parameter values along with the 
SDSS + data from JS]. The figures are shown for the case Q m = 0.27, h = 0.71, but 
the qualitative features of the spectrum are independent of the exact value of Q m . The 
spectra are normalized to match the matter power spectrum of the ACDM concordance 
model at k = 0.01 /iMpc" 1 . 

From the figures it is clear that the matter power spectrum is drastically modified, 
much like in the Chaplygin gas scenario studied in |13j . Even a small deviation from 
the ACDM values (n = 0, q = 1) leads to either strong oscillations or growth in the 
spectrum, depending on the signs of n and q — 1. This is explained by the non-zero 
sound speed that seeds the non-standard evolution of 5 my k due to the cj.k 2 / a 2 -term in 
Eq. (|TK|) . We know generally that perturbations in a comsological fluid with c 2 > 
will fluctuate for scales smaller than the Jeans length, which is the scale when pressure 
start to dominate over gravity. While for fluids with c 2 < 0, they will be unstable and 
grow exponentially for the same scales. A more thorough discussion on perturbations 
for fluids with non-zero sound speed can be found in 

Expanding the sound speed for the MPC model (}2*2*|) in the small quantities n and 
q — 1, we find that they determine the sign of c 2 . For n > and q — 1 < c 2 will be 
negative and hence give rise to exponentially growing perturbations, while for n < 
and q — 1 > c 2 will be positive and the perturbations will oscillate. This is exactly the 

+ http://www.sdss.org 
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Figure 4. The matter power spectrum in the MPC model without baryons. for small 
deviation of q (a) and n (b) from the ACDM values of q = 1, n = 0. Also shown are 
data points and error bars from the SDSS project. 



behaviour we see in Fig. This behaviour for large k can also be understood directly 
from Eq. ()18|) . In the large k limit we can ignore all but the k 2 term in the factor 
that multiplies the zeroth derivative of the density perturbation. Then the differential 
equation simplifies to 

c 2 k\ , x 

(23) 



<5 m , k + if (2 - 3c 2 s ) <5 m , k + -^-6 : 



m,k 







We can get a feeling for how the density perturbation behaves in this limit by looking 
at the case where the time variation of a and c~ is relatively slow. Then both a and c 2 
will be approximately constant, and the first derivative term in the differential equation 
will vanish. The differential equation reduces then to the common harmonic oscillator 
equation, with the following solution 

<W(*o) ~ e~ icMo - (24) 
From this we see that a positive c 2 will give oscillating perturbations, while a negative 
c 2 gives exponentially growing perturbations. We expect the perturbations to exhibit 
such behaviour for scales smaller than the Jeans length. The Jeans length is given by 
the expression 

(25) 



For the Cardassian model the Jeans length today can be written as 



X 2 

A J,0 



]ir 2 n 



i n 



c 



3 h 2 i + n c 

The pressure starts to dominate when 
k > = kj 



in 



l)q 



1 + nQc 



A 



(26) 



(27) 
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Let us look explicitly at the Cardassian model with parameters n = and q = 1.00001. 
Using the expression for the Jeans length in Eq. ([26)1 . we find that kj « 0.6h/Mpc. 
According to our simplified analysis, the perturbations should then start to grow 
exponentially when k > 0.6h/Mpc. This agrees fairly well with what we see in Fig. 

Hi. 

In jlHj it was argued that the very tight constraints in the Chaplygin case arising 
from the strong deviations in the matter power spectrum could be relaxed by considering 
the effect of baryons. Even though the constraints are relaxed slightly, the overall effect is 
very strong and the fluctuating Chaplygin gas scenario remains tightly constrained |13j . 
even without considering the overall normalization of the spectrum, which tightens the 
constraints even more. Here, we expect a similar situation to apply but for completeness 
consider the effect of adding the baryons to the system. 

For k > H the evolution of linear density perturbations of the Cardassian and 
the baryon fluids are given approximately by a coupled set of second order differential 
equations. Taking the equation of state parameter and the adiabatic sound speed of 
matter to be zero, these differential equations can be written as ([El EH): 

h + H{2 - 3(2u - c 2 s ))5 k - (Ih 2 (7w - 3w 2 - 6c 2 s ) - (^) 2 )^ 

- 3 -H 2 (l + w) p6k + Pb6b - k =0 (28) 
2 P + Pb 

and 

4 + 2i/4-^±^ = 0, (29) 

2 p + Pb 

where 5 is again the density perturbation of the Cardassian fluid, 5b the density 
perturbation of the baryon fluid and w and c 2 are the equation of state parameter 
and adiabatic sound speed of the Cardassian fluid. 

We solve the equations as before and the resulting matter power spectrum is shown 
in Fig. 03 We see that adding the baryons to the system does improve the situation as 
the oscillatory behaviour disappears, and the main effect is a reduction in small-scale 
power. The constraints are still very strong, but bearing in mind that we have neglected 
other effects which can result in less power on small scales, e.g. neutrino masses and a 
scale-dependent bias between matter and galaxy clustering, this simple analysis cannot 
rule out these models completely. However, as shown in the integrated Sachs- Wolfe 
effect in the CMB leads to very strong constraints on models with c 2 > 0. In the cases 
where c 2 < 0, the exponential blow-up is still present, and is strongly disfavoured by 
the galaxy power spectrum. 

Here, considering the effect on erg, does not tighten the constraints like in the 
Chaplygin gas case. This can be seen by comparing how the linear growth is changed 
in the two scenarios: the Chaplygin gas case is shown in [46J and the MPC case in [24J. 
Linear growth in the Chaplygin gas scenario changes radically for very small deviations 
from the parameter values corresponding to the ACDM model, whereas for the MPC 
model this is not the case. We have checked that this also holds for the parameter values 
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Figure 5. The matter power spectrum in the MPC model with baryons for small 
deviation of q (a) and n (b) from the ACDM values of q = 1, n = 0. The oscillatory 
behaviour has now vanished for n < and q > 1 but still very small deviations from 
the ACDM model worsen the fit to SDSS data. 

considered here by calculating the linear growth in the MPC scenario for n < 0, q < 1. 
Linear growth is modified, but not as radically as in a Chaplygin gas dominated universe 
and the LSS constraints are much tighter. Therefore adding baryons gives somewhat 
more freedom in choosing the parameters in the MPC scenario, but one should bear 
in mind the strong constraints from the ISW effect [H]. As is clear from Fig. the 
parameter space is highly constrained and the model offers no good alternative to the 
ACDM model. 

5. Conclusions 

We have considered the observational constraints arising from current CMB, SNIa and 
LSS data on commonly considered models within the Cardassian framework, the original 
Cardassian model and the Modified Polytropic Cardassian model. Both models were 
introduced as an alternative explanation to the acceleration of the universe. 

Considering that the models only have an effect on the background evolution, CMB 
and SNIa observations constrain the allowed parameter space. We find that, within 
the approximation that only background evolution is modified, using the full CMB 
spectrum constrains the original Cardassian model significantly and that the model offers 
no advantage over the ACDM model. Adding a new degree of freedom by extending 
the consideration to the Modified Polytropic Cardassian model relaxes the constraints 
somewhat but still the parameter space is more constrained than what one finds using 
the supernovae data only. The constraints indicate that a somewhat non-standard model 
is slightly preferred over the ACDM model but not significantly. 
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Viewing the Cardassian models as arising from a interacting fluid constrains the 
scenario critically. Considering the growth of fluctuations along with the SDSS data, we 
find that any deviation from the ACDM model is strongly disfavoured, with or without 
including baryons in the calculation. The fluid interpretation of the Cardassian model 
is hence effectively excluded as a viable alternative to the ACDM model. This is further 
supported by recent work [H], where the CMB TT power spectrum in the Cardassian 
model with fluctuations is demonstrated to be modified strongly on large scales. 

In short, even small deviations from the ACDM model within the fluid 
interpretation of the Cardassian model are strongly disfavoured by current data. If 
on the other hand, only background evolution is modified, the data indicates a slight 
preference over the ACDM model but the difference is not significant. 
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